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HYDROGEN-OXYGEN CHEMICAL REACTION KINETICS 
IN ROCKET ENGINE COMBUSTION 
by Martin Hersch 
Lewis Research Center 

SUMMARY 

Hydrogen-oxygen reaction times and concentration histories of chemical species 
during reaction were calculated for rocket combustor conditions. Calculations were 
made for oxidant-fuel weight ratios of 1 and 10 and initial reactant temperatures of 
1200° to 2500° K at a chamber pressure of 20 atmospheres (2.0x10® N/sq m). The 
reaction time varied from about 0. 01 second at 1200° K to a few microseconds at 2500° K. 
Calculations were made by using a numerical integration program and an analytical 
solution. The reaction mechanism used included five chain branching reactions and 
three recombination-type reactions. 


INTRODUCTION 

Various steps in the rocket combustion process are atomization, vaporization, gas 
phase mixing, and chemical reaction. The slower steps control combustion efficiency. 
Some steps may also drive combustion instability, if their completion time approximates 
the wave time. Therefore, it is important to know the completion time for the various 
steps, the overall combustion' time, and instability wave times. Combustion times for a 
typical rocket engine are reported in reference 1. High-frequency instability times are 
treated in reference 2. 

Vaporization is often the slowest and, therefore, rate controlling step. Accordingly, 
it has been well studied (ref. 3). Gas phase mixing may also be a slow step and is con- 
sidered in reference 4. Less information is available for chemical reaction rates in 
rocket combustion. Often, it is simply assumed that the reactions take place infinitely 
fast in comparison with other steps and thus may be neglected. This assumption may not 
always be valid. For example, if the slower steps, such as vaporization and mixing, are 
made highly efficient, then it would be important to know the chemical reaction times. 
Furthermore, combustion instabilities are known to cause large temperature and pres- 
sure gradients. Since chemical reaction rates are sensitive to temperature and pressure, 
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the assumption of infinitely fast reaction rates under these conditions may also be 
invalid. 

Calculations of hydrogen -oxygen (H-O) kinetics for reaction limited ram -jet and 
rocket combustor conditions are presented in references 5 to 7. In these studies the 


following reaction mechanism was used: 

H + 0 2 sr O + OH I 

O + H 2 r H + OH H 

H 2 + M = 2H + M HI 

h + oh+m: h 2 o + m iv 

0 2 +Mz20+M V 

H 2 + OHzH+ h 2 o VI 


where M is any third body. 

In references 5 to 8 the initial temperature was varied from approximately 1000° to 
3000° K at pressures of 1 to 65 atmospheres (1x10® to 6. 5x10® N/sq m). Brokaw, using 
an analytical solution for reactions I, H, VI, VH, and VHI (ref. 9), shows that the 
following reactions are also important, particularly at low temperatures: 

ho 2 + h 2 s h 2 o 2 + h vn 

h + o 2 + m =ho 2 + m vm 

The importance of these two reactions increases rapidly as the initial temperature de- 
creases. It would therefore be desirable to consider a mechanism which includes all 
eight reactions. 

Calculations for this system of eight reactions are presented in this study. Numer- 
ical and analytical techniques are used. The effects of the initiation reaction 
Hg + Og — 2 OH is also considered with the use of the analytical technique. 

The numerical method considers reversible reactions, concentration changes of all 
species, and temperature change during reaction. The analytical method is restricted 
to constant H 2 and 0 2 concentrations, nonreversible reactions, and constant temperature. 
Reaction times are presented as functions of initial temperature and oxidant-fuel weight 
ratios O/F at a constant pressure. Histories of atom and radial concentrations during 
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the reaction are also shown. Calculations were made for initial temperatures ranging 
from 1200° to 2500° K and O/F values of 1 and 10 at a chamber pressure of 20 atmos- 
pheres (2.0X10® N/sq m). 


CALCULATIONS 
Numerical Solution 

The numerical integration program used for this study is that developed under NASA 
contract and described in reference 8. The assumptions made in this program are 

(1) The mixture is one of thermally perfect gases. 

(2) Flow is inviscid throughout. 

(3) Transport properties can be neglected. 

(4) Internal energies are in thermal equilibrium. 

(5) The law of mass action applies throughout. 

Backward rates are calculated, with microscopic reversibility assumed, from the 
ratio of forward rate constant to equilibrium constant. Equations of mass, energy, and 
state are solved simultaneously with the differential equations for specie concentration 
changes due to chemical reactions. 

Input data for the numerical program are temperature, pressure, species concentra- 
tions, forward rate constants, and equilibrium constants. Parameters for the forward 
rate constants are given in table I, and parameters for the equilibrium constants in 
table n. The equilibrium constants were calculated from the data of reference 10 and 
similar unpublished Lewis data. 

The reactions, including those proposed by Brokaw (ref. 9), used for the numerical 
integration program are 


H + 0 2 ;0 + OH 

In 

O + H 2 rr H + OH 

Hn 

2H + M =t H 2 + M 

mn 

+ OH + M2 H 2 0 + M 

IV n 

20+M:0 2 + M 

Vn 

H 2 + OH = H + H 2 0 

VIn 
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TABLE I. - FORWARD REACTION RATE PARAMETERS IN EXPRESSION 51 


K = DT n exp (-E A /RT) 


Reaction 

number 

D, 

3 

cm , mole, sec 

n 

Activation 

energy, 

E A’ 

cal/mole 

Refer- 

ence 

Reaction 

I 

5.64X10 13 

0 

15 100 

5 

H + o 2 rr O + OH 

n 

1.2X10 13 

0 

9 200 

5 

O + H 2 5= H + OH 

in 

2.0X10 18 

-1.0 

0 

7 

2H + M - H 2 + M 

IV 

3. OxlO 19 

-1.0 

0 

11 

H + OH + M = H 2 0 + M 

V 

1.85X10 17 

5 

0 

5 

20 + M r 0 2 +M 

VI 

4.2X10 12 

. 5 

10 000 

5 

H 2 + OH = H + H 2 0 

vn 

5.4X10 11 

0 

24 000 

9 

HOg + H 2 sr H + H 2 0 2 

b vm 

c ^Fh 2+ 1 .7,5F 02 ) 8 .® <1014 

0 

-1 280 

12 

H + 0 2 + M =r H0 2 + M 

IX 

2. 5 OxlO 12 

0 

39 000 

13 

H 2 + 0 2 =r 2 OH 


^Where K is rate constant, D is preexponential constant, T is temperature in °K, n is 


temperature exponent, E A is activation energy, and R is universal gas constant. 
^Third-body efficiencies obtained from unpublished Lewis data. 
c Where F is mole fraction. 

TABLE H. - EQUILIBRIUM CONSTANT PARAMETERS IN EXPRESSION 3 
K C =D C T c exp(F c /T) 


Temperature range, 500° to 4600° K.j 


Reaction 

number 

D c> 

cm, mole, sec 

E c 

F c> 

°K 

Reaction 

I 

4.0748X10 2 

-0.4103 

-8. 6625X10 3 

h+o 2 :0 + oh 

n 

1. 8764 

.02212 

-9. 1499X10 2 

o+h 2 :h + oh 

m 

.27658 

-6.444X10' 3 

5.2421X10 4 

2H + M 2 H 2 + M 

IV 

5.9727X10" 3 

.2779 

6.0344X10 4 

H + OH+Mr H z O + M 

V 

1.2722X10' 3 

.4262 

6.0169X10 4 

20+M^0 2 + M 

VI 

2. 1570xl0~ 2 

.2845 

7. 9234X10 3 

h 2 + oh:h + h 2 o 

vn 

.6614 

-.0464 

-7. 6370X10 3 

ho 2 + h 2 :h + h 2 o 2 

vm 

6.9973X10 -2 

.2799 

2.3858X10 4 

h + o 2 + m: ho 2 + M 


^here K c is equilibrium constant, and D c , E c , and F c are equilibrium 
constant coefficients. 



vnn 


H0 2 + H 2 sr h + h 2 o 2 

H + 0 2 + M 3T H0 2 + M V TTTn 


Analytical Solution 


An analytical solution is obtained with isothermal conditions and constant reactant 
(H 2 and 0 2 ) concentration for the system assumed: 


H + O,- 


K, 


o + h 2 - 


K<- 


K 0 


H 2 + M 


K c 


0<; + M • 

Ca 


K c 


H 2 + OH 


K r 


ho 2 + h 2 . 


K c 


H + 0 2 + M 


H0 2 + M 


K_8 


-O + OH 
H + OH 

► 2H + M 
-2 O + M 

► H + H 2 0 
- H + H 2 0 2 

HOg + M 
H + 0 2 + M 


la 

Ha 

IHa 

Va 

Via 

VHa 

vma 

-vnia 


Reactions TTTa and Va are written in the direction of dissociation because during initiation 
and induction radical formation is of concern. Reaction IVn is omitted from the analyti- 
cal solution because it results in a term containing products of radicals, which results in 
a nonlinear equation. This reaction, however, is of little importance during the induction 
period. Reaction VHI is considered in both directions in the analytic solution. 

The differential equations expressing the growth rate of intermediate species are 
(for the case of reaction VH! in the forward direction) 


* 
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(la) 


= K 1 C 0 2 C H - K 2 C H 2 C 0 + 2K 5 C 0 2 C M 


dC 


H 


dt = ' KlC °2 C H + K 2 C H 2 C 0 + 2K 3 C H 2 C M 


+ K 7 C H2 C H02 - K 8 C 02 C m C h + K 6 C H2 C OH 


dC 


ho 9 

— — = - K 7 C H2 C HC>2 + KgC 02 C M C H 


~ K 1 C 0 2 C H + K 2 C H 2 C 0 ‘ K 6 C H 2 C OH 

where the C is concentration, K is a reaction rate constant, and t is time. 
The general solution for these differential equations is (ref. 9) 


A X t 

C 0 =£ A i e + a o 
i=l 


* A.t 

C H = £ B i e 1 + b H 
i=l 


X.t 


C H0 o = E C i e 1 + 


'HOr 


i=l 


X.t 


C 0H = £ D i el+ dOH 

i=l 


where A^, B., CL, and are constants with dimensions of concentration, and 


.-1 


constant with dimensions of second . 

The terms a Q , bjj, c r q , and 

2 

given in reference 9. For example, 


The terms a Q , b H , c r q , and d OH , which result from initiation reactions, 

2 


(lb) 

(lc) 

(l d) 

(2a) 

(2b) 

(2c) 

(2d) 

X is a 
are 
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( 3 ) 


a o = 


K 3 C H 2 C M + K 5 C 0 2 C M 

V7 


Standard matrix techniques may then be used to determine the values of X. Then, with 
the boundary condition that the radical concentrations are zero at time zero, the pre- 
exponential constants may be calculated. It then remains to calculate the growth of 
H 2 0 and H 2 <D 2 . The differential equations for these are 


dC 


H 9 0 

dT" = K 6 C H 2 C OH 


(4a) 


dC 


H 2°2 

dt 


K 7 C H 2 C H0 2 


(4b) 


Since equations (4) contain terms which are now known, they may be integrated directly, 
with the boundary conditions of zero concentration at time zero. The detailed analytic 
solution was prepared and programmed for a high-speed digital computer by David A. 
Bittker of the Lewis Research Center. 


RESULTS AND DISCUSSION 

Results are presented for an oxidant-fuel weight ratio O/F of 1 and 10, a chamber 

e 

pressure of 20 atmospheres (2. 0x10° N/sq m), and initial temperatures ranging from 
1200° to 2500° K. 

The principal reason for studying chemical kinetics in rocket combustion is to 
determine the reaction time. This time is the time required to convert vaporized and 
mixed propellants to burned combustion products. For the hydrogen-oxygen system the 
reaction takes place in three distinct stages. The first stage is a very short initiation 
period during which small concentrations of atoms and radicals, such as H, O, and OH, 
are produced. This period is followed by the induction period, characterized by constant 
exponential growth of atoms and radicals, under nearly isothermal conditions. During 
induction a plot of the logarithm of concentration against time is a straight line. 

Very late in the induction period atom and radical concentrations become large 
enough to cause self -heating due to three -body recombination reactions. The temperature 
then greatly increases, atom and radical growths depart from constant exponential 
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increase, and the system attains equilibrium. This final period may be referred to as 
the postinduction period. 

The numerical and analytical methods are somewhat complementary. The numerical 
calculations cannot be started at the true beginning of the reaction, with the assumption 
that the reaction starts with zero atom and radical concentrations. On the other hand, 
the analytical program can be started for the conditions of zero atom and radical con- 
centrations at zero time. The analytic method, however, gives no information which 
can be used to define the end of induction. 

The OH concentration plotted against time during initiation and induction for various 
initial temperatures is shown in figure 1. Temperature change is also shown. The 
solid and dashed curves represent numerical and analytical results, respectively. The 
lower left ends of the curves show the very short period of initial formation of atoms 
and radicals. The curves then show the relatively long induction period characterized 
by constant exponential growth under isothermal conditions. The numerical results show 
that the induction period, as defined by constant exponential growth, continues, in some 



Figure 1. - OH growth and temperature rise. 
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Figure 1. - Concluded. 

cases, well into the region of temperature increase. In this study the induction period 
is considered ended when the atom and radical concentrations deviate from constant 
exponential growth. This deviation can be determined, in this study, by the numerical 
method only. Although the analytical assumptions are violated when the system is 
nonisothermal, the limit to which the analytical results can be continued is the point at 
which the concentration curves deviate from a straight line, as shown in figure 1. This 
point has therefore been selected as the end of induction and is plotted at zero time. 

Numerical computation times increased drastically with decreasing temperature. 
Therefore, at 1200° K numerical results were obtained only near the end of induction 
and thus are not shown in figure 1. 

The OH concentration at the start and end of induction are shown in figure 2. The 
concentration at the start of induction was determined by the intersection of the ex- 
trapolated slopes of the initiation and induction portions of the curves, such as shown 
in figure 1(a). The induction period is often considered ended when the OH concentration 

Q 

reaches the threshold of experimental detectibility, 10 gram-mole per cubic centi- 
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1200 1400 1600 1800 2000 2500 3000 

Initial temperature, °K 

Figure 2. - OH concentration at start and end of induction. 



Figure 3. - Temperature rise during postinduction period (numerical 
solution). 



meter (ref. 9). This criterion appears to be suitable only for an initial temperature 
near 1500° K. The induction time may now be determined by using these limits of OH 
concentration shown in figure 2. 

Since in this study the total reaction time, rather than only induction time, is im- 
portant, the postinduction time must also be investigated. Some idea of the post- 
induction time may be obtained from the curves of figure 3, which shows the temperature 
rise for various initial temperatures as a function of time. These results show that for 
the temperature range of about 1500° to 2500° K the postinduction time is on the order 
of a few microseconds. Though only partial results are shown for 1400° K, apparently 
at this temperature, the postinduction time may be on the order of perhaps 30 micro- 
seconds. 



Time 

Solution 

Oxidant-fuel 

ratio, 

O/F 

- Induction 

Numerical 

1 

- Induction 

Analytical 

1 

•- Induction 

Analytical 

10 

- Postinduction 

Numerical 

1 


(approximate) 


1200 1400 1600 1800 2000 2500 

Temperature, °K 

Figure 4. - Effect of initial temperature on induction and postinduction time. 


The induction and approximate postinduction times are shown in figure 4. Numerical 
and analytical methods show nearly identical results for the induction time. For an 
initial temperature higher than about 1700° K, the induction time is small in comparison 
with the postinduction time and therefore may be neglected. At lower temperatures the 
postinduction time, though, becomes increasingly small in comparison with the induction 
time. The induction time increases rapidly with decreasing temperature for temper- 
atures below 1600° K. In this region the induction time increases exponentially with the 
reciprocal of initial temperature. 

The effect of O/F on induction time was investigated by determining the OH con- 
centration at the end of induction for an O/F of 10 by using the numerical program. 
Calculations at an O/F of 10 were made for temperatures between 1500° and 2500° K. 
Detailed calculations at this O/F for the entire induction period or for lower temper- 
atures were not made because of the long computational times. The OH concentration 
at the end of induction for an O/F of 10 is shown in figure 2. At high temperatures the 
OH at the end of induction at an O/F of 10 is about 10 times that at an O/F of 1. At 
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temperatures lower than 1700° K the OH concentration at the end of induction is approx- 
imately the same for both O/F conditions. 

The induction time based on these numerically calculated OH concentrations at the 
end of induction for an O/F of 10 is shown in figure 4. The curve was extended to 
1200° K by assuming that the OH concentration at the end of induction for an O/F of 10 
was the same as for an O/F of 1. These results indicate that at higher temperatures 
the induction time decreases when the O/F is increased. At lower temperatures the 
effect of O/F on induction time appears to be small. 



ca 



Figure 5. - Relative product species during induction period (normal- 
ized to OH concentration). 


Some idea of the composition during induction is given by figure 5. During initia- 
tion relative concentrations of the atom and radical species shift with time. During 
induction, however, the relative concentrations remain nearly constant with time. 
Relative concentrations of O, H, HC^, HgO, and HgC^ are shown, relative to OH 
concentration, during induction as a function of temperature. Almost perfect agreement 
is obtained between the numerical and analytical results. The relative concentrations 
of O and H appear nearly constant during induction over most of the temperature range. 
Concentrations of HC> 2 , HgO, and HgOg are much greater than OH at low temperatures. 
At higher temperatures the concentrations of O, H, HOg, and HgO are one or two orders 
of magnitude greater than OH. At 2500° K the I^Og concentration during induction is 
nearly two orders of magnitude less than that of OH. 

All numerical results presented thus far were obtained by starting the computations 
with initial concentrations of atoms and radicals determined analytically, because the 
numerical program cannot be started with zero atom and radical concentrations. The 
effect of starting numerical calculations with arbitrary atom and radical concentrations 
is shown in figure 6. The dashed lines are concentration histories when the starting 
concentrations were calculated analytically. The solid curves are concentration 
histories for arbitrary initial concentrations. These results show that if the initial 
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Figure 6. - Comparison of concentration histories for arbitrary with 
analytically determined initial concentrations. Temperature, 

2000° K; pressure, 20 atmospheres (2. OxlO 6 N/sq m) ; oxidant- 
fuel ratio, 1. 


concentration of a species is too large, the concentration will remain essentially constant 
until the correct relative concentration is reached. It will then increase as it would 
have, had its initial concentration been correct. 

It is possible to include additional initiation reactions in the mechanism. The most 
important of these reactions is perhaps the reaction: 


H 2 + 0 2 - 2 OH 

The effect of this reaction on the time to reach a given OH concentration is shown in 
figure 7. Including this reaction decreases the reaction time by about two-thirds. At 
higher temperatures, where the postinduction time is controlling, this shortening of the 
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Figure 7. - Effect of reaction H 2 + 0 2 ^ 20H 
on time to attain arbitrary OH concentra- 
tion (analytical calculations). Oxidant- 
fuel ratio, 1. 


induction time would have no effect on the overall reaction time. This reaction might be 
of more importance at lower temperatures. It can be shown, however, that for rocket 
combustion studies this is not a significant effect. 


SIGNIFICANCE OF CHEMICAL REACTION 
TIMES IN ROCKET COMBUSTION 

Liquid-propellant rocket-engine combustion involves many steps, atomization, 
vaporization, gas phase mixing, and chemical reaction. The importance of any steps 
in the combustion process may be related to its speed or completion time. The slower 
steps will be more important than the rapid ones. The combustion dead time, an in- 
dication of the time required for the overall process, is reported in reference 2 for a 
typical H 2 -0 2 rocket engine. It varied from 0.25 to 2.25 milliseconds for a wide range of 
operating conditions. High-frequency instability must also be considered. Reference 1 
shows that wave times for high-frequency instability range from about 0. 1 to 1 milli- 
second. Thus, the time range of interest in liquid rocket combustion ranges from 0. 1 
to several milliseconds. Therefore, if chemical reaction times become greater than 
about 0. 1 millisecond, they must be considered in studying rocket combustion. 

This study shows that chemical reaction time increases very rapidly with decreas- 
ing reactant temperatures if the temperature is below about 1700° K. At 1600 K the 
reaction time is about 10 microseconds, which may be considered infinitely fast in com- 
parison with other steps in rocket combustion. At about 1400° to 1500° K the reaction 
time increases to about 10"^ second and thus enters the region of interest in rocket 
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combustion. Above 1700° K the time decreases to several microseconds, and is fairly 
insensitive to temperature. 

It is beyond the scope of this study to investigate the heating of unreacted propellants 
in the combustion chamber. Nor can it be determined here what the gas temperature is 
at the start of reaction. Presumably, though, the propellants in actual rocket combus- 
tors are substantially heated soon after injection by convection, radiation, and recircu- 
lation. The results of this study indicate that if the unreacted propellants are rapidly 
heated to about 1400° K, then chemical kinetics can be neglected as a rate controlling 
mechanism in steady-state combustion for a chamber pressure of 20 atmospheres 
(2.0x10® N/sq m). 

Combustion instabilities, however, are known to cause large temperature and pres- 
sure gradients. If the instability wave reduces the temperature to below about 1500° K, 
then reaction times would greatly increase and become approximately equal to the wave 
period. Under these conditions the reaction kinetics could become an instability driving 
force. 


SUMMARY OF RESULTS 

Hydrogen-oxygen reaction times and concentration histories were calculated for 
rocket combustor conditions. Calculations were made for oxidant-fuel weight ratios 
of 1 and 10 and initial reactant temperatures ranging from 1200° to 2500° K at a 
chamber pressure of 20 atmospheres (2.0x10 N/sq m). A numerical integration pro- 
gram and a programmed analytical solution were used. The following principal results 
were obtained: 

1. Hydrogen-oxygen reaction times varied inversely with initial reactant temper- 
ature, ranging from a few microseconds at 2500° K to approximately 0. 01 second at 
1200° K. Since processes requiring times less than milliseconds can be neglected in 
determining performance efficiencies or combustion characteristics this study indicates 
that chemical reactions can be neglected if the injected propellants are heated above 
approximately 1400° K. 

2. When the initial propellant temperature was higher than about 1700° K, the re- 
action time was primarily the postinduction or heat-release period. At lower temper- 
atures the reaction time was the isothermal induction time. 

3. Induction time decreased with increasing oxidant -fuel mixture ratio. This effect 
of the oxidant-fuel mixture ratio on induction time increased with increasing temper- 
ature . 
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4. The analytical method presented permitted very rapid and accurate calculations 
of induction times and concentration histories for the hydrogen-oxygen system. How- 
ever, in order to calculate induction times additional information must be available to 
define the end of the induction period. The numerical results presented provided this 
information. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 1, 1967, 

128-31-06-03-22. 
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